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Introduction 

Research In fundamental studies in geodynamics has continued in a 
number of fields including seismic observations and analysis, synthesis 
of geochemical data, theoretical investigation of geold anomalies, 
extensive numerical experiments in a number of geodynamlcal contexts, 
and a new field — seismic volcanology. 

Summaries of work in progress or completed during this report 
period are given below. Abstracts of publications submitted from work 
in progress during this report period are attached as an appendix. 

Analysis of Long-Period Surface Waves Excited 
by Mt. St. Helens Eruption on May 18, 1980 
— Terrestrial Monopole? 

The eruption of Mt. St. Helene on May 18, 1980 excited 

long-period seismic surface waves which were recorded clearly by the 
IDA, SRO and ASRO stations. Some examples are shown in Figure la, b and 
c. In general. Love waves are larger than Rayleigh waves. The spectra 
of these surface waves are computed and equalized to the source. Figure 
2 shows the amplitude spectra of Love and Rayleigh waves as a function 
of azimuth around the source. Both Love and Rayleigh waves show a 
two-lobed radiation pattern with a nodal direction in about N10°E and 
E10°S for Love and Rayleigh waves respectively. 

These radiation patterns have been Interpreted by Prof. Kanamorl 
and Jeff Given using normal mode theory. It can be shown that no 
double-couple source can produce the radiation patterns shown in Figure 
2. The radiation patterns for a single horizontal force, a single 
vertical force and an isotropic source are computed and are 
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schematically shown In Figure 3. Comparison of Figure 3 with Figure 2 
Immediately suggests that an almost horizontal single force In either 
NIO^E or S10*^^ direction Is a good representation of the source. 
Kanamorl and Given performed a detailed Inversion of the source spectra 
at periods of 100, 120, 150, 200 and 256 sec to determine the geometry 
of the source. On the basis of this inversion they conclude that the 
force Is nearly horizontal (with E15° from the horizontal) and Its 
horizontal projection is In the direction S5®W. 

By using the source spectrum and the synthetic seismograms they 
determined the time history of the source. They assumed the source time 
function in the form: 

s(t) “ 1/2(1 -coS'£ii)) 0 ^t-2x 

- 0 t>2t 

and computed synthetic seismograms for various values of t (Figure 4). 

Comparison of the synthetics with the observed records Indicates that t 

« 75 sec Is appropriate. From the amplitude spectrum, the peak value of 

1 ft 

the force Is determined to be 10 dyne. 

This result provides Important constraints on the mechanics of 
volcanic eruption. They plan to present this result at the coming Fall 
AGU meeting in San Francisco, and a full paper which describes the 
details of the analysis is now In preparation. 
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Upper mantle anisotropy 

There is a growing body of evidence that ouch of the mantle is 
anisotropic with respect to the propagation of seismic waves* Professor 
Anderson, working with Professor Dziewonski of Harvard and with students 
at the Seismological Laboratory, is investigating this phenomenon. 

The presence of anisotropy is not only of theoretical interest. 
If the upper mantle is, in fact, anisotropic then isotropic inversion of 
seismic data will result in erroneous structures because of improper 
parameterization. Such important geodynamlcal problems as the presence 
and nature of a low-velocity zone and the depth extent of differences 
between oceans and continents depend critically on the validity of the 
assumption of Isotropy. 

Since intrinsic anisotropy requires both anisotropic crystals and 
preferred orientation, the anisotropy of the mantle contains information 
both about the mineralogy and the flow pattern in the mantle. For 
example, olivine, the most abundant upper mantle mineral, is extremely 
anisotropic for both P-wave and S-wave propagation. It apparently is 
easily oriented by the ambient stress or flow field. Olivine-rich 
outcrops show a consistent preferred orientation over large areas. In 
general, the seismically fast axes are in the plane of the flow with the 
a-axis, the fastest direction, pointing in the direction of flow. The 
b-axis, the minimum velocity direction, is generally normal to the flow 
plane. On a global basis, Anderson and Dziewonski find that the average 
anisotropy in the upper 200 km of the mantle is that which would result 
from the horizontal shearing associated with pl&te motion. There is 
evidence that under ridges, the sense of anisotropy is reversed. 


consistent with upwellins flow there. 

Dynamic calculations of geold enoaalles 

Geoid anomalies ere usually interpreted in terms of density 
contrasts buried in a rigid Earth. Theoretical work by graduate student 
Mark Richards and Professor Hager shows that for a more realistic 
viscous earth, density contrasts cause flow and stresses which deform 
the top surface and any internal boundaries, as illustrated in Figure 5. 
This boundary deform&tlo. aas au important effect on the total geoid 
anomaly - in fact, for a uniform viscosity earth, a negative geoid 
anomaly results from a positive density contrast. The total anomaly can 
be an order of magnitude larger than the anomaly from the density 
contrast alone, as shown in Figure 6. This is because the stresses 
causing deformation of the surface fall off less rapidly with depth than 
the potential from the density contrast. 

These calculations have been generalized to include an earth with 
varying viscosity and with chemical layering. The magnitude and sign of 
the total anomaly is a strong function of the properties of the layers. 
Finite element code for geodynamics applications 

The rheology of the mantle is complicated, as is the geometry in 
many situations of interest in geodynamics. Numerical experiments with 
realistic geometries and rheologies are ess ntial to test hypotheses and 
to give insight into the physics of geologic processes. The finite 
element technique is best suited for application to gsodynamics 
problems. 

Professor Hager, along with Dr. Arthur Raefsky, Visiting 
Associate, has spent a substantial portion of this period developing a 
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specialized finite element code tailored for solution of problems of 
Interest in geodynamics. This code, except for cosmetic revisions, is 
complete and a manual for its use is being written. Much of the 
numerical work reported below was done using this code. 

Ridge Processes 

Two very important observations are that oceanic crustal thickness 
is independent of spreading rate and that the entire crustal section is 
formed completely essentially at the ridge axis. The first observation 
has been Interpreted in terms of a passive model of a ridge, in which 
upwelllng under ridges fills in the void left by spreading plates. In 
this model, isothermal upwelllng beneath ridges is a passive response 
to, rather than a driving mechanism for plate motions. 

Professor Hager has calculated the flow pattern beneath ridges 
according to the assumptions of the passive model and using a realistic 
nonNewtonian rheology. Streamlines are shown in Figure 7. These 
streamlines are more widely spaced beneath ridges than they are beneath 
plates, indicating that upwelllng occurs over a wide region. Melt will 
rise vertically, so this passive model is inconsistent with the 
requirement that the crust be formed at the ridge axis. 

Professor Anderson arrived at the conclusion that ridges are 
active, not passive, from geochemical considerations. In his model, 
dlaplrs of pier it Jc melt rise from an olivine-eclogite layer at a depth 
of 220 km. 

It is essential in Anderson's model that the diaplrs rise rapidly 
enough that crystals not settle. He argued previously that the effects 
of nonNewtonian rheology will lead to rapid ascent. Hager and Raefsky 
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have computed more detailed finite element modela to quantify the 
vifect. Aacent velocltlea of meters/day are to be expected. 

The flow pattern within a rlalng dlaplr will tend to keep crystals 
from settling. This flow, as well as the ascent velocity of a dlaplr, 
depends on Its shape. Graduate student Cathryn Allen has determined 
using numerical experiments that dlapirs evolve Into spheres. She Is In 
the process of determining crystal settling criteria for rising dlapirs. 
Earthquake migration and rheology of the asthenosphere 

Four large earthquakes occurred In a seismic gap of category 2 
(McCann et al. , 1979) near the Loyalty Is. in the New Hebrides. At 

3;25 on October 24, 1980, an event with Mg ■ 6.7 initiated the sequence. 
Three events, llg ■ 6.7, 7.2, and 6.5 followed this event on the next 
day. 

Graduate student John Vldale and Professor Kanamorl investigated 

this sequence by using the PDE data, f lrst*nnotlon data from WWSSN 

records, and long-period surface waves recorded by IDA, SRO, and ASRO 

stations. The flrst<nBotlon data constrained one of the nodal planes. 

With this constraint, inversion of Rayleigh and Love-wave spectra at 256 

seconds determined the other nodal plane. The mechanisms of all four 

events were almost pure thrust on a plane dipping 20° and striking 

parallel to the local strike of the New Hebrides trench. The first-day 

aftershocks indicate an initial rupture zone of about 2,000 km which Is 

27 

consistent with the escimated seismic moment of 3 x 10 dyne-cm. This 
cluster of activity Is shown in Figure 8. During the next two days, the 
aftershock activity expanded to an area of 110,000 to 20,000 km^ in the 
directions both along and perpendicular to the trench. The expanded 
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zone of seismicity is slso shown in Figure 8. This pattern suggests 
that the initial rupture zone represented a zone of increased strength 
(i.e. asperity), and the stress change due to failure of this asperity 
subsequently migrated outward. The events on the %rest side of the 
trench, the largest of which appears to have been a normal event from 
its first motions, were almost undoubtably triggered by the thrust 
events. Two similar thrust sequences in the New Hebrides have caused 
similar aftershock normal events. 

We are using this data set, one of the cleanest of its kind, to 
investigate the dyanamlc response of the lithosphere to stress loads 
from large thrust events. We expect to be able to place tighter 
constraints on the coupling between the lithosphere and asthenosphere 
and on the rheology of the asthenosphere. This data should provide a 
useful test of the rheological model previously developed under this 
grant by Anderson and Minster, 

Seismic coupling at subductlon zones 

Graduate student Larry Ruff and Professor Kanamorl have previously 
shown that M^^, the maximum characteristic size of earthquakes at 
subductlon zones, is well correlated with convergence rate and plate age 
(Figure 9), in an attempt to understand the physics governing this 
correlation, they have examined earthquake source parameters for 
representative events in the northern Pacific, 

They find that the maximum earthquake size is governed by the 
asperity distribution on the fault plane. Large events are 
characterized by rupture of large asperities, while smaller events are 
characterized by rupture of multiple smaller asperities (see Figure 10), 
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It it llktly chat asperity site it related to the state of stress 
in the slab, which is in turn related to slab trajectory. Graduate 
student Jia Jung Zhang is studying the relation between slab trajectory 
and stress in slabs using finite eleaent oodels constrained by observed 
gravity, topography, and slab trajectory. 

The global tsantle flow pattern oust influence slab trajectory. 
Professor Hager has computed pressures from the global return flow for a 
variety of models. He finds that for models with flow extending 
throughout the mantle, maximum earthquake size correlates well with this 
pressure field. 

Sel8m.,city and composition of volcanic rocks in island arcs 

Plate subduction is generally considered to be the ultimate cause 
of large earthquakes and volcanism at Island arcs. The activity of 
large earthquakes can be related to various parameters of plate 
subduction. Here we measure seismic activity by the magnitude of the 
largest characteristic earthquake for e&ch subduction zone. 

There is a general correlation between and plate convergence 
rate V, and and the age of the subducting plates, T, as shown in 
Figure 9. This correlation Is not unexpected because of reflects the 
strength of inter-plate coupling, and V and 1/T are considered to be 
proportional to the strength of plate coupling (Ruff and Kanamorl, 
1980). On the other hand, various parameters for arc volcanic rocks 
have been associated with other characteristic features of island arcs. 
Among them are the ratio CA/TH (Cal-Alkali Series/Tholeiltic Series) 
(Miyashiro, 1974), Silica Index Q (Suglmura, 1968), Alkali (K 2 O, Na20 
etc.) (Kuno, 1966; Dickinson and Hatherton, 1967; Dickinson, 1968; 
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Nielson end Stolber, 1973; Suglsakl 1976; U1 and Araaaki, 1978) and 
H 2 O conte'<ts (Sakuyama, 1981). These parameters have been correlated 
with th>' depth of the Uadati-Benlof f sons; crustal thickness, Bouguer 
gravity anomalies and the plate convergence ra'.e. 

One such example is shown in Figure 11 (Coulon and Thorpe, 1981). 
These correlations have been interpreted in terms of: (1) different 
depth of fractional crystallization, (2) different degree of partial 
iselting, (3) different degree of contwalnation by crustal rocks during 
magma ascent. 

In these mechanisms, important elements are the density contrast 
dp between the crust and the magma and the viscosity n in the upper 
mantle both of which control the depth of partial melting, fractionation 
and the residence time of the magma in the upper mantle and the crust. 

Since represents the strength of coupling, it may be 
proportional to the horizontal compressive stress La in the upper plate. 
Although La is of the order of 10 to 100 bars, it is comparable to the 
deviatoric stress responsible for ascent of mantle diapirs. Numerical 
models of diapirs show o:hat La affects the velocity of magma ascent, 
and works in the same manner as Ap and u in controlling the composition 
of volcanic rocks at island arcs. 

In order to explore tisis possibility, we examined the correlation 
between and K 2 O (at Si02 37 to 1 Z) content, and and the 
percentage of basalt for various island arcs. The results are shown in 
Figures 12 and 13. Despite fairly large scatter in the data, 
correlations seem to indicate that the regional variation of the 
horlzonal co^>ressive stress represented by is at least partially 
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rcsp'tnsiblc for th« roglonal variation of volcanic rocks. 

In the course of this study, we found a rather serious difficulty 
in obtaining a uniform data set for volcanic rocks. This study is still 
at a preliminary stage, and we are collecting and compiling more data 
from various sources. 

Flow and stress in subduct ion tone regions 

Subduction zones provide excellent places for comparison of 
numerlcil models and observations. Seismicity in Benioff zones reveals 
the location of subducted slabs and reflects the state of stress within 
slabs. 

Earthquakes In subduction zones are not uniformly distributed. 
The variation of earthquake activity with depth is presumably related to 
1) variations of viscosity with depth, 2) variations of density with 
deptn or 3) some time dependent process that involves plate ages, 
convergence rates, temperature and kinetics. The various parameters of 
each subduction zone, e.g. , number of earthquakes, energy release, age 
of lithosphere, convergence rate, total length of subduction zone, dip, 
age of associated arc, presence of back-arc spreading, are being 
assembled. 

The number of earthquakes and the energy release as a function of 
depth are particularly interesting. (tost Benioff zones exhibit an 
exponential die off of activity from the surface to about 250 km. When 
old lithosphere, > ''5 m.y., is subducted there is generally a secondary 
peak in activity between 300 and 670km. There are clearly at least two 
populations of events. In some cases there is a more gradual decay of 


activity in the upper mantle followed by an abrupt drop at 200-250 km 
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depth. 

Graduate atudent Marius Vaasiliou and Professor Anderson have 
classified different regions into three categories, depending upon their 
selsnicity (see Table and Figure U). In Class I, seismicity ceases 
near 250km, with more or less linear decay to that depth. In Class 2, 
the deepest seismicity is near 550 km. This clacs is divided into two 
subclasses. In the first (2a), the decay is roughly linear above ~250 
km below which there is a gap, and a separate peak centered around 
4400-500 km. In the second (2b), the decay is also roughly linear above 
~250 km. (although less so than in 2a), after which there ia a break in 
slope, and a slower decay to ~550 km. Class 3, where the deepest 
seismicity is arou. t 670 km, is likewise divided into two subclasses. 
In (3a), there is detached deep peak, with shallow seismicity extending 
to ~250 km. In (3b), there is continuous seismicity, with a break in 
Che decay slope at ~250 km, a low level of seismicity below that, and a 
resurgence at >00-600 km. Regions in Class 1 tend to involve subduction 
nf fairly young crust, and those in Claases 2 and 3 tend to involve 
subduction of older crust. Each class may represent a stage in the 
evolution of subduction. Plausible sequences of evolution might be 1 > 
2a 3a or 1 ♦ 2b ♦ 3b. 

The seismicity presumably reflects stress in the slab. Vassili'Aj 
and Hager are computing finite element models of flow and stress near 
subducted slabs. An example of stress in a vertically sinking slab is 
shown in Figure 15. There is a roughly linear decrease in acress in the 
upper part of the slab, which could explain the aeismicity obaervations 
nicely if the number of sarthquakea depended exponentially on the stress 
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(c.g., ttogl, 1962). Th< Incrctit In strctt in th« low«r pare of tha 
•lab la tuggaatlva of Class 3 raglons. Modals ara balng axcandad to 
invastigsts ths affset of slab dip and coupling tha slab Co cha surface 
Plata. 

Profassor Hagar has invascigacad cha daforoacion of cha 670 ka 
saisaic discontinuity in chasa aodals. convaction is confiiMd to cha 
uppar asncla by chanici«l stratif icacion, slabs should panacrata 25-1^0 
ka into tha lowar asncla below cha average depth of this disconcinuicy. 
Saisaic daca collactad by Research Fallow Ichiro Nakanishi should help 
Co resolve this iaportsnt question. 

Tha location of Banioff sona salsaicicy conscrslns c‘ s asncla flow 
patcarn since it reveals the position of subduccad slabs. Profassor 
Hager axtandad his previous global flow aodels to include the affects of 
trench migration and increased strength of subducted slabe. Ha found 
that considering these affects does not alter his previous conclusion 
that for Che Newtonian mantle viscosities dk.rivad from postglacial 
rebound studies, flow confined to the upper rsntlc does not aatch the 
observed dips of subducted slabs. Hsntle-wida circulation passes this 
isqiortant test. 

Saisaic array studies 

Studies dene by gradusta student Marianne Ualck using the 
Caltech-U.S.G.S. SCARLET array are providing data which proaisa to 
provide useful constraints on our understanding of geodynanical 
processes. In a study with Bernard Hinstar, Ualck studied lateral 
heterogeneities in P-wave velocities beneath Southern California. She 
was sble to explain P-wave travel time anoaalies in teras of a chin lens 
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model at a depth of 150 km. One possible Interpretation of this anomaly 
is that it is due to alignment of olivine crystals in a downwelling flow 
beneath the Big Bend in the San Andreas fault. Hager is developing the 
mathematical formulation of a three-dimensional flow model to test this 
hypothesis. 

Ualck is now using the array in a study with Professor Anderson to 
study upper mantle structure beneath western Mexico and Gulf of 
California. The velocity gradient between the 400 km and 650 km 
discontinuities is well constrained and similar to that for northern 
Eurasia. 

An Important, unanswered question in Southern California 
geodynamics is the mechanism for the support of the Transverse Ranges, 
which reach elevations in excess of 10,000 ft. The extent to which 
these mountains are supported by crustal thickness variations is being 
explored by graduate student Thomas Hearn in conjunction with Minster. 
Hearn is using delay times of arrivals, shown in Figure 16, to 
examine crustal thickness variations. The presence of a crustal root 
beneath the Sierra Nevada is confirmed. Thin crust, less than 20 km, 
occurs in the Salton Trough and extends into Arizona. Thickness of over 
30 kis are prevalent throughout the Transverse Ranges. 
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Distribution of Regions in the 3 Classes 
(see text) 
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Figure 1: 


Figure 2: 


Figure 3: 


Figure 4; 


Figure 5: 


Figure 6: 


Figure 7: 


Figure Captions 

Records of long period surface waves excited by the eruption 
of Mt. St. Helens, May 18, 1980. 

Anplltude spectra of Love (G) and Rayleigh (R) waves from the 
Mt. St. Helens eruption as a function of azimuth. 

Radiation patterns for a single horizontal force, a single 

vertical force, and an Isotropic source. 

Synthetic seismograms for various source time histories 
compared with that observed for the Mt. St. Helens 
eruption. 

Cartoon representing two of the three contributions to geold 
anomalies in a fluid earth. A dencity contrast a within the 
mantle Induces flow. This flow results in stresses at the 
surface and core-mantle boundary, which deform. This 
deformation also causes a geold anomaly. A third effect, 
added stresses resulting from the gravity field from a 
acting on the mantle, is not shown. 

Ratio of fi'^oT ® function of depth to the density 

contrast and spherical harmonic degree for a uniform 
viscosity mantle. The total potential Is of opposite sign 
than that for a rigid earth, 6V^. For large I the magnitude 
of the total geold anomaly Is much greater than that which 
would result in a rigid earth. 

Flow vectors and streamlines for a passive ridge model with 
nonMewtonlan (power law three) rheology. The various dashed 
lines depict the volume of the mantle sampled for a 6 km 
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thick basaltic crust generated by partially melting a given 
fraction of the mantle. 

Figure 8: Migration of seismic activity after the New Hebrides 

earthquake of October 24, 1980 shown on a bathymetric map of 
the area. The heavy line cHitlines the zone of initial 
activity. An outward migration of activity is evident. The 
cluster of events in the fourth frame (10/26) are normal 
fault events seaward of the trench. 

Figure 9: The correlation between earthquake size and: convergence 

rate and oceanic lithosphere age. The number at each 
subduction zone is the characteristic My (or Mg) value, and 
the regression plane is shown as contours of constant My. 
The broken line encloses the subduction zones (in the lower 
left corner) where there is either confirmed or suspected 
back-arc spreading. The important trend is that the largest 
earthquakes occur in subduction zones with a faster rate and 
younger age (see Ruff and Kanamorl, 1980, for details). 

Figure 10: The difference in the source time histories and inferred 

asperity distributions of the 1964 Alaskan and 1963 Kurile 
Islands earthquakes. The seismograms are shown at top as 
the solid traces. The dashed traces are the synthetic 
seismograms corresponding to the deconvolved source time 
functions shown beneath the seismograms. While the Kurile 
Islands time function is composed of a multiple event 
sequence, the Alaskan time function represents a smooth 
rupture spreading over a large area. Notice the large 
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Figure 11: 


Figure 12: 


Figure 13: 


Figure 14: 


Figure 15: 


Figure 16: 


difference In the seismic moments also* As discussed In 
Ruff and Kanamorl (1981), the time functions of these events 
can be Interepreted as the spatial dlstrilHition of 
asperities. The epicenters (stars) and fault areas are 
shown below, and the Inferred asperities are Indicated as 
the hachured areas. The Alaskan earthquake occurred In a 
region of huge asperity length scale. 

Correlation between crustal thickness and crustal parameters 
(from Coulon and Thorpe, 1981). 

Correlation betwen maximum earthquake size and potassium 
content. 

Correlation between maximum earthquake size and basaltic 
component of volcanlcs. 

Logarithm of number of earthquakes versus depth, with a depth 
Interval of 20 km., for earthquakes occurring since 1964 

with m^ > 4 for some regions to lllutrate the different 

classes referred to In the text, (a) Ryukyus, an example 
of Class 1. (b) ■ Japan, Class 2a. (c) « Java, Class 2b. 

(d) > Tonga, Class 3a. (e) ■ Peru, Class 3b. 

Stress magnitude and velocity In slab for a finite 

element model of a vertical sinking slab. 

Contours of delay times (In units of 0.1 sec) relative to 

station ABL. One second of relative delay time results from 

a relative thickness variation of 10 km, given uniform 
crustal velocity. Negative delay times Indicate a thinner 


or slower crust 
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A GLOBAL GEOCHEMICAL MODEL FOR THE EVOLUTION OF THE MANTLE 

Don L. Anderson 

Seiemological Laboratory, California Institute of Technology 
Pasadena, California 91125 


Abstract . Basalt, eclogite, and 
harzburgite, differentiation products of the 
Earth, appear to be trapped in the upper mantle 
above the 670 km seismic discontinuity. It is 
proposed that the upper mantle transition region, 
220 to 670 km, is composed of eclogite, or 
olivine eclogite, which has baen derived from 
primitive mantle by about 20Z partial melting and 
that this is the source and sink of oceanic 
lithosphere. The remainder of the upper mantle 
is garnet peridotite, or pyrolite, the source of 
continental basalts and hotspot magmas. This 
region is enriched in Incompatible elements by 
partial melts or hydrous and CO 2 rich metasomatic 
fluids which have depleted the underlying layers 
in the L.I.L. elements and L.R.E.E. The eclogite 
layer is internally heated. It may control the 
convection pattern in the upper mantle. Material 
can only escape from this layer by melting. The 
insulating effect of thick continental 
lithosphere may lead to partial melting in both 
the peridotite and eclogite layers. Hotspots 
and ridges would then mark the former locations of 
continents. Most of the basaltic or pyroxenltlc 
fraction of the oceanic lithosphere returns to 
the eclogite layer. 

In this model plate tectonics is intermittent. 
At a depth of 150-220 km the continental thermal 
anomaly triggers kimberlite and carbonatite 
activity, alkali and flood bcsalt volcanism, 
vertical tectonics and continental breakup. Hot- 
spots remain active after the continents leave, 
building the oceanic Islands. Mantle plumes rise 
from a depth of about 220 km. Midocean ridge 
basalts originate in the depleted layer below this 
depth. Material from this layer may also be 
displaced by subducted oceanic lithosphere to fom 
back-arc basin magsms. 


Introduction 

Although convection plays an important role in 
plate tectonics and heat transport in the Earth, 
it has not succeeded in homogenizing the mantle. 
Magmas are still being produced from mantle 
reservoirs which have remained separate for the 
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Abstract . The high surface tenperature of 
Venus laplles a penMnently buoyant lithosphere 
and a thick basaltic crust. Terrestrial style 
taetcmles with deep subductlon and crustal re- 
cycling is not possible. Orerthlckmed basaltic 
crust partially aslts Instead of converting to 
ecloglte. Because aentle nagnaa do not have 
convenient access to the surface the ^"Ar 
abundance In the ataosphere should be low. 

Venus nay provide an aiudog to Archean tectonics 
on the Barth. 


Introduction 

The surface of Venus la about 450 K wrmmr 
than the surface of the Earth. This affects the 
buoyancy, themal expansion, therval conduct Id ty 
and, hence, the themal evolution and ultimate 
fats of the lithosphere. The pressure in the 
upper mantle of Venus Is at least 12X less than 
at equivalent depths In the Earth's upper mantle. 
Thus, the depth of melting, the locations of 
upper mantle phase changes and the viscosity of 
the upper mantle will be different for the two 
planets. The buoyancy and thermal properties of 
the llthos^ere control the style of plate 
tectonics and the associated time and length 
scales. Most discussions of the comparative 
tectonics between Earth and Venus address only 
the differences in viscosity. 

The oceanic lithosphere of the Barth cools as 
It ages and It eventudly becomes denser than the 
underlying mantle. This Instability develops 
after about 40 m.y. (Oxburgh and Paramentler, 
1977). Once cooling roaches a depth of 50 km 
the garnet pyroxenlte or ecloglte tnnsformatloos 
may make a substantial contribution to the 
negative buoyancy of the llthos^ere. Continents 
and oceanic plateaus resist subductlon because 
of their thick low-density crusts. As fate tf 
the terrestrial lithosphere, therefore, depends 
on both chemistry and temperature. Ae purpose 
of this paper Is to Investigate the Implications 
for Vmus teetonlds of Its hl^ surfsce 
temperature. He will conclude that the surface 
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Su— ry 

Itocrople Earth aedala art uaabla to prorida unifora fits to 
tha groaa Earth noraal soda dataaat or, In nany caaaa, to raglonal Lora 
and laylai|h waya data. Anisotropic inrorsion provides a good fit to 
tha data and indicates that tha upper 200 ka of tha aantlt is 
anisotropic. Tha nature and aagnituda of tha required anisotropy, 
•oraovar, is sinilar to that found in body wave studies and in studies 
of ultranafic saaples froa the upper aantle. Pronounced upper Mntle 
low-velocity sones are characteristic of aodels resulting froa isotropic 
inversion of global or regional datasets. Anisotropic aodels have aore 
nearly constant velocities in the upper aantle. 

Noraal aode partial derivatives are calculated for a transversely 
Isotropic Earth aodel with a radial axis of syaaetry. For this type of 
anisotropy there are five elastic constants. The two shear type aoduli 
can be deterained froa the toroidal aodes. Spheroidal and Rayleigh 
aodes are te::sitive to all five elastic constants but are aainly 
controlled by the two coapressional-type aoduli, one of the shear-type 
aoduli and the reaaining, aixed-aode, aodulus. The lack of sensitivity 
of Rayleigh waves to coapressional wave velocitisr Is a characteristic 
only of the isotropic case. The partial derivitives of the horizontal 
and vertical coaponents of the coapressional velocity are nearly equtl 
and opposite in the region of the aantle where the shear velocity 
sensitivity is the greatest. The net co^ressional wave partial 
derivative, at depth, is therefore very saall for isotropic 
perturbation. Coapressional wave anisotropy, however, has a significant 
effect on Rayleigh wave dispersion. Once it has been established that 
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transverse anisotropy Is Important it is necessary to invert for all 
five elastic constants. If the asiauthal effect has not been averaged 
out 1 !»ore general anisotropy aay have to be allowed for. 

Introduction 

There is a growing body of evident ^ that ouch of the upper auintle 
auy be sotropic to the propagation of seisaic waves. The early 
evidence was the discrepancy between Rayleigh wave and Love wave 
dispersion (Anderson, 1961, 1966, Anderson and Rarkrider, 1962) and the 
asiauthal dependence of oceanic velocities (Raitt et al., 1971, Hess, 
1964). Asiauthal variations have now been documented for isany areas of 
the world (Blbee and Shor, 1976, Baaford, 1977) and the Raylelgh-Love 
discrepancy is also widespread (McEvilly, 1964, Aki, 1968, Forsyth, 
1975, Yu and Mitchell, 1979, Schlue and Knopoff, 1977). Shear wave 
birefringence, a aanifestation of anisotropy, has also been reported 
(Anderson, 1966, Him, 1977, Masataks et al., 1980). The degree of 
anisotropy varies but is typically about 5Z. It is not known to what 
depth the anisotropy extends but soise of the data requires it to depths 
as great as 200 ka. 

Anderson (1966) proposed that the "discrepancy" between aantle 
Rayleigh and Love waves could be explained if the vertical P and S 
velocities in the upper asntle were 7-8Z less than the horisontal 
velocities. He also showed that aodels without an upper aantle 
low-velocity xone, such as the Jeffreys aodel, could satisfy the 
dispersion data if the upper aantle was anisotropic. The Love-Rayleigh 
discrepancy has survived to the present and aodels have been proposed 
which have SV in the upper aantle less than SH by about 4Z. These more 
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THE PHYSICS OF CREEP AND ATTENUATION IN THE MANTLE 


Don L Andenon and J. Bernard Minster 

Seisnological Laboratory, California Institute of Technology, Pasadena, California 91 125 


Abstract. Dislocations contribute to both seismic wave attenuation 
and steady-state creep in the mantle. The two phenomena involve quite 
different strains and characteristic times but they can both be understood 
with simple dislocation models. The most satisfactory model for creep 
involves the glide of dislocations across subgrains and rate limited by the 
climb of jogged dislocations in the walls of a polygonized network. The 
jog formation energy contributes to the apparent activation energy for 
creep, nuking it substantially larger than the activation energy for self- 
diffusion. The theory leads to either a or o* law for creep rate, 
depending on the length of the dislocations relative to the equilibrium 
spacing of thermal jogs. 

Attenuation in the mantle at seismic frequencies is probably caused by 
the glide of dislocations in the subgrains. Kink and impurity drag can 
both contribute to the glide time constant. The kink-formation, or Peierls 
barrier, model for dislocation gUde appears to be a low-temperature, high- 
frequency mechanism most appropriate tor pure systems. A small amount 
of impurity drag brings the dislocation glide characteristic time into the 
seismic band at upper-mantle temperatures. 

The attenuation and creep behavior of the mantle are related through 
the dislocation structure. Discussion of the various possible mechanisms 
is facilitated by casting them and the geophysical da*a in terms of a pre- 
exponential characteristic time and an activation energy. The relaxation 
strength is an additional parameter that can be used to identify the 
attenuation mechanism. Mobile dislocations in subgrains, rather than cell 
walls, have the appropriate characteristics to explain the damping of 
seismic waves in the upper mantle. The grain boundary peak may be 
responsible for attenuation in the lithosphere. 

1. Introduction 

The puipose of this paper is to investigate the role of dislocations in 
both aeep and attenuation. The high temperature creep of most 
crystalline solids is controUed by dislocation climb or glide plus climb 
(Weertmann. 1^)68]. For a wide variety of metallic, ionic, and polar 
crystals the activation energies for aeep and self-diffusion are esentially 
identical, indicating that the motions of dislocations are rate-limited by 
self-diHusion. Furthermore, in the stress range of about 10"^ to 10"* G, 
where G is the shear modulus, the stress exponent of the creep rate is 
about 3, which is consistent with dislocation climb or the glide of jogged 
dislocations. Both of these mechanisms of steady-state creep require self- 
xxxxxxxxxxxxxx 
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Abstract 

Subducted slabs penetrate a convectlng, actively flowing mantle. 
Some of the variations observed among subduction tones may result from 
differences in the global flow pattern and thereby provide constraints 
on mantle flow models. We compute global flow models by solving for 
flow in a viscous sphere subject to the observed plate velocities as 
boundary conditions. Density contrasts assumed appropriate for the 
surface plates and subducted slabs are included. Models are considered 
in which flow extends throughout the mantle or in which flow is confined 
to the upper mantle. 

For both types of model, all subduction zones which are currently 
undergoing back-arc spreading are characterized by a shear in the 
asthenosphere tending to pull the overriding plate toward the trench. 
All island arcs which are not undergoing back-arc spreading, with the 
possible exception wi Sunda, have the opposite sense of shear 
contributing to back-arc compression. Instantaneous flow velocity 
vectors match observed Benloff zone dips for the model which allows 

mantle wide flow, but not for the upper mantle model considered, which 

22 

has a viscosity of 10 p. The maximum earthquake magnltuda observed 
for island arcs is well correlated with age of the subducted slab and 
pressure gradient between the trench and back-arc region if flow extends 
throughout the mantle. The effect on particle trajectories of trench 
migration is not important if subducted slabs sink faster than trenches 
migrate. Variation in slab dip and seismic coupling can be explained by 
global flow if flow extends throughout the mantle. If flow Is confined 
to the upper mantle, its effective viscosity must be substantially less 
than lO^^p. 
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Abstact . We present a correction to a paper by 
Okal and Anderson (1975) about multiple ScS 
travel-time anomalies. We have reanalyzed data 
for ScS 2 surface bounces In the South Atlantic 
Ocean. From these data an ScS 2 ~S residual of 23.6 
seconds was found by Okal and Anderson (1975). 
This corresponded to an SCS 2 surface bounce point 
under Trlndade Island and was Inferred to be due 
to very slow upper mantle associated with the 
Trlndade hot spot. The analysis we present here 
Invalidates this conclusion. The nature of the 
upper mantle under Trlndade Is an open Issue. 

Introduction 

It Is possible to use ScS 2 ~ScS time anomalies 
to estimate S travel-time heterogeneities In the 
upper mantle be’ow the surface bounce point of 
SCS 2 . The time difference between SCS 2 and ScS or 
S Is obtained by cross-correlation of the two 
signals. The residual Is the difference from the 
predicted time using Jef f reys-Bullen (JB) tables 
(1958). This method was used by Okal and Anderson 
(1975) -hereafter referred to as O&A- with data 
from the WWSS long period network. O&A obtained 
travel time residuals for bounces under continents 
and oceans. Among their numerous data, one Is of 
special Interest both because It gives the largest 
delay and because It corresponds to the unique 
case of a surface bounce under a 'hot spot'. 
Indeed O&A found an 11.8 seconds one-way 
travel-time residual corresponding to an 8082 
surface bounce point under Trlndade Island In the 
South Atlantic Ocean (see figure 1). Trlndade's 
volcanlsm Is as young as 1 million years ,and 
extremely silica undersaturated, suggesting origin 
from a hydrous mantle (Oversby , 1971 ; Valenclo and 
Mendla, 1974) . This led O&A to propose that the 
large 8 delays which they observed were due to 
reflection In a very slow upper mantle, possibly 
Indicating a higher degree of partial melting 
under the hot spot. In this paper we present a 
detailed analysis of the data O&A used to reach 
this conclusion. We show that for the South 
Atlantic data points O&A used the SV component 
Instead of the SH component. In fact, the stations 
used are on a node of the SH radiation pattern. 
Besides the fact that waveform distortion might 
occur on the SV component, we show that In this 
particular case the first arrival corresponds to 
the Interference of three phases; SKS,S and ScS. 
Moreover we examine seismograms for stations 
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Abstract 

Seismic coupling has been used as a qualitative measure of the 
"Interaction" between the two plates at subductlon zones. Kanamorl 
(1971) introduced seismic coupling after noting that the characteristic 
size of earthquakes varies systematically for the northern Pacific 
subductlon zones. A quantitative global comparison of many subductlon 
zones reveals a strong correlation of earthquake size with two other 
variables: age of the subducting lithosphere and convergence rate. The 
largest earthquakes occur in zones with young lithosphere and f>’st 
convergence rates, while zones with old lithosphere and slow rates are 
relatively aseismic. Results from a study of the rupture process of 
three great earthquakes indicate that maximum earthquake size is 
directly related to the asperity distribution on the fault plane ( 
asperities are strong regions that resist the motion between the two 
plates). The zones with the largest earthquakes have very large 
asperities, while the zones with smaller earthquakes have small 
scattered asperities. This observation can be translated into a simple 
model of seismic coupling, where the horizontal compressive stress 
between the two plates is proportional to the ratio of the summed 
asperity area to the total area of the contact surface. While the 
variation in asperity size is used Uo establish a connection between 
earthquake size and tectonic stress, it also implies that plate age and 
rate affect the asperity distribution. Plate age and rate can control 
asperity distribution directly by use of the horizontal compressive 
stress associated with the "preferred trajectory" ( l.e. the vertical 
and horizontal velocities of subducting slabs are determined by the 
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plate age and convergence velocity). Indirect Influences are many, 
including; oceanic plate topography and the amount of subducted 
sediments. 

All of the subduction zones are apparently uncoupled below a depth 
of 40 km, and wc propose that the basalt to eclogite phase change in the 
down-going oceanic crust may be largely responsible. This phase change 
should start at a depth of 30-35 km, and could at least partially 
uncouple the plates by superplastic deformation throughout the oceanic 
crust during the phase change. 
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ABSTRACT 


Spatial averaging of short-period teleseismic P-wave arrivals 
across SCARLET constrains upper mantle lateral velocity variations 
beneath southern California. We compare dT/dA, back-azimuth and 
averaged arrival time estimates determined from the entire network for 
154 events to the same parameters derived from small subsets of SCARLET. 
The resulting relative array diagrams and plots of "net subarray delays" 
computed for IVi overlapping subarrays exhibit large anomalies that vary 
smoothly as a function of subarray position, 1‘^plying substantial 
lateral heterogeneity under the network. Patterns of mlslocation 
vectors for more than 100 sub^rrays permits identification of an 
East-West striking high-velocity antiform beneath the Transverse Ranges. 
Further constraints on both the lateral configuration and depth of the 
anomaly are obtained through thin lens analysis of the net delay data. 
The best-fit lens — at 150 km depth — features low velocities beneath 
the Salton Trough as well as a concentrated high-velocity zone between 
116.3° - 118. 3°W and 33.7° - 34.8°N, consistent with the relative array 
diagrams end the PKP-delay contours of Hadley and Kanamori (1977). 
Therefore the bulk of the high-velocity material observed by Hadley and 
Kanamori may in fact lie at 150 km rather than 40 - 100 km. This 
possibility prompts renewed speculation regarding the relationship of 
this body to the plate tectonic history of the region. 



MODELS OF FLOW AND STRESS BENEATH RIDGES 


B. H. Hager (Selsmo Lab, Caltech, 

Pasadena, California 91123) 

Two classes of models of flow beneath ridges 
are presented. Kinematic (passive) models of 
flow driven by the motion of the plates are 
calculated analytically for both Newtonian and 
power law rheologies. Plots of the deviatoric 
stress tensor are shown. For both rheologies 
the deviatoric stress below the ridge is zero, 
leading to the conclusion that shear heating 
from kinematic flow beneath ridges is of 
negligible importance in the generation of 
partial melts. Vertical flow for n=3 creep 
takes the form of a broad plug with the maxi- 
mum upwelllng velocity a factor of two smaller 
than that for Newtonian creep. Power law 
creep leads to a shallower sampling of the 
mantle beneath ridges, tending to Isolate the 
deeper mantle. 

In the .second type model, we investigate the 
buoyant upwelling of a partially molten diapir. 
A model with Newtonian rheology is used as a 
basis of comoarlson for a series of models 
which employ more realistic rheologies. The 
effects on rheology of stresses generated by 
expansion due to Increased melting are Included. 
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DEFORMATION OF SEISMIC DISCONTINUITIES 
AND THE SCALE OF MANTLE CONVECTION 

Bradford H . H.? gcr 

Arthur Rarfsky (Both at: Sei sinological 

Laboratory, Caltech, Pasadena, CA 91125) 

The stresses from the horizontal density con- 
trasts in a converting system cause defornwition 
of the system boundaries. Readily observable 
examples <if this effect are the ~A km elevation 
of the seafloor at midore.iiilc ridges and the ~6 
km depression at treiiclios. Numerical models 
show that for convection confined to the upper 
mantle .•Imi ,i ifie entire weight of the subducted 
llthosplien- is compensated by deformation of the 
surface and fi70 km discontinuity. 

For a given stress, tlie deformation of an in- 
terface is inversely projior t ional to the density 
contrast .u ross it. If tlie 67t) km discontinuity 
were a clninical discontinuity, not a pliase boun- 
dary, the deformation for a given stress would 
be 5-10 tint'.: that for the same stress at the 
surface. 'Ihi distribution of this deformation 
depends upon the rheology of the lithosi'here and 
mantle. Sl.ib strength and decoupling by the 
asthenosphere cause stresses at 670 km to he 
much large.' than those supported by deformation 
of the upper surface. Depression of the 670 km 
discontinuity bineatb slabs would be at least 25 
km and would likely reach 150 km. A flat dis- 
continuity Is compatible with a phase change 
with negligible Clapeyron slope but not with a 
chemical discontinuity; in particular, cessa- 
tion of seismicity at the global average depth 
of the 670 ku di scont Ircjity is not compatible 
with a layered nvintlc. 

Ph.isc changes, coupled wltii deform.it ion of 
chcLiical discontinuities, can lead to finite 
amplitude instabilities in mantle convection 
wfiicb result in episodic overturns on a large 
and small scali-. Tills provides a meclianism for 
mixing between layers. Examples Include exten- 
sive mc‘llirig caused by pulling a chemical layer 
above Its solidus and massive overturn of a 
layered mantle If a subducted slab penetrates 
deeply enough to undergo x pliase change. 
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MODELS OF FLOW AND STRESS NEAR SUBDUCTED SLABS 


B. H. IlaRer (Selsmological Laboratory, 252-21, California Institute of 

Technology, Pasadena, California 91125) 

A. Racfsky (Jet Propulsion Laboratory, Pasadena, California 91109) 

R. J. O'Connell (Dept, of Geological Sciences, Harvard University, 

Cambridge, Massachusetts 02138) 

One important observation constraining the dynamics of subduction is 
that plates are bent well beyond their elastic limits upon subduction, yet 
are then restraighteaed, rather than remaining curled. Another fundamental 
observation is that slabs do not in general have dips close to vertical. 
These observations imply that the net effect of tractions acting on the 
slab from viscous flow within the mantle is comparable to or larger than 
the negative buoyancy of the slab and much larger than the slab’s 
resistance to bending. The angle of subduction, as well as the state of 
stress within the slab, is strongly influenced by the mantle flow. 
Observations of slab dip and state of stress, in combination with models 
of mantle flow, can help to constrain mantle rheology and provide a 
better understanding of the dynamics of plate motion and subduction. 

Flow in the mantle is influenced by both the large-scale global flow 
associated with plate motions and small-scale flow driven by local density 
variations. Previously calculated kinematic models of the global flow 
driven by the (imposed) motion of the surface plates show that for models 
which allow return flow deep in the mantle, the instantaneous flow 
vectors match the observed dips of deep Benioff zones (Hager and O'Connell, 
Tectoiiophysics , 19 3). These models arc highly simplified, however, in 
particular, neglecting the increased density of the slabs, the supposed 
increased strength of slabs relative to the asthenosphere, and the 
temporal evolution of plate boundaries. 

New global models have been calculated which include the density 
contrasts appropriate to the subducted slabs, located on the basis of 
observed seismicity. For models which allow flow deep in the mantle, the 
match between instantaneous velocity vectors and Benioff zone dip-; is 
Improved; velocities beneath subduction zones are close to the velocities 
of the surface plates, so presumably flow vectors closely match slab 
trajectories in magnitudes as well as dips. For the model in which flow 
Is confincJ to the upper mantle, the constraints of return flow dominate 
and there is little change in velocity from the previous kinematic models. 
Flow cuts across the deep Benioff zones in the Sunda, Tonga, and Middle 
America regions. 

Consideration of reasonable trench migration velocities has little 
effect on the models which allow deep mantle flow, although the shallow 
flow models are changed substantially. Although it is conceivably possible 
to manipulate shallow-flow models to match the observed dips of Benioff 
zones, the fact that the simple deep flow nodels match the dips so vrell 
is evidence that the flow associated with plate motions extends deep into 
the mantle. 
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llie effects of t)ie likely lelatlvely hlplu-r strenp.tli of the slab are 
inves;tit’,nt.o«.l ur.lnp, a scries of two-d iiTiensional finite clemont models. 
Plausible contrasts in the effective viscosity of the slab relative to 
the asthenosphere have surprisinp.ly little effect on the orientation of 
the flow lines. The local state of stress is more sensitive than 
velocities to variation in the rheolc>p.y of the earth's inter lor. 

Deviations from the Ruff and Kanamori (PEri, 1980) empirical relation 
bet ween seismic coupling;, convergence velocity, and plate age appear to 
be related to regions where the global return flow is complicated; seismic 
rones where coup! iirg is anomalously weak tend to be regions where, the 
global flow exerts pressure oir the upper surface of subducted slabs. Deep 
mantle flow may effect seismic coupling at shallow depths as well as slab 
dip and stress. 


TONGA 


CATHUS RIGID 



Section.'; through flow models iir the Tonga region. The models on the left 
have a 10^^ p lower mantle, while the models on the right have an 
effectively rigid lower mantle. In all llu* model:;, the obser vi'd jrlate 
velocities are imposed as boundary conditions. The lower row of models 
has a density contrast appropriate to a subducted slab in the reg,ion of 
earthquake hypocenters (marked by boxes). The d i sjil acement s j>lotted 
represent the Inst antaneou:: place velocities ext rai'olnt ed for 7 My. 



TIME-TERM MAPPING OF THE MOHO DISCONTINUITY IN 
SOUTHERN CALIFORNIA 

Thomas M. Hearn ( Selsmological Laboratory, 

California Institute of Technology, Pasadena, 

California 91125) 

J. Bernard Minster (Systems, Science and 

Software, P.O. 1620, La Jolla, CA 92038) 

The time-terra method for the analysis of 
refraction data (Willmore and Bancroft, 1960) 
has been applied to Pn arrivals at the southern 
California seismic array. A combination of more 
than 30 earthquakes and 175 stations produce the 
data set of over 2000 Pn arrivals. The 
Inversion finds a mean Pn velocity of 8.0 km/sec 
for the entire data set, however variations of 
.2 kra/sec from this value are known to exist in 
southern California. The delay times computed 
are equlvaleat ,„to the function 

/q^( (l-V(z) '/Vji )^'^/V(z))dz where V(z) is the 
crustal velocity, is the Pn velocity, h Is 

the Moho depth and z is depth. The variation of 
over 1.5 seconds in the delay times is 
indicative of about 15 km of relief in the 
crust-mantle boundary. 

A thinning of the crust to less than 20 kin 
occurs in the Salton Trough with the thin crust 
apparently extending Into southwest Arizona. A 
previous refraction study in Arizona (Warren, 
1969) indicates a crustal thickness of 20 km 
southwest of Gila Bend. This suggests that the 
crustal thinning associated with the Gulf of 
California begins in the southeast California - 
southwest Arizona region. Crustal thickness in 
the Channel Islands averages about 24 km while 
thicknesses of over 30 km are prevalent 
throughout the Transverse Ranges. This is 
consistent with the higher Bouguer gravity 
anomalies of the Channel Islands. The high 
delay times of the Los Angeles and Ventura areas 
are due to tfie thick sedimentary sequences of 
the associated basins. The presence of a 
crustal root beneath the Sierra Nevada is 
confirmed by events at Mammoth which demonstrate 
a delay of more than .5 seconds for paths 
travelling directly along the bathollth. 
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MKCIIANISM 01' MT. ST. IIKLENS ERUPTION DETERMINED 
FROM EONG-PERIOD SURFACE WAVES 


Ui roo K 

Jcffroy W. Given (both at: Selsnulogiral l.abora- 
tory, California Institute of Tccbnolocy, 
Pasadena, California 91125) 

The eruption of Mt . St. Helens on May 18, 
1980 excited lonj;-period (100 to 260 see) surface 
waves, and high quality seisinop.ram.s were recorded 
at many IDA, SRO and ASRO stations. Tne.se are 
prob.ably the first dat.i of this kind obtained for 
a voleiiiie event, and provide a unicjue opportu- 
nity I" study tlie inecbanism of eruption. The 
amplitude i.idiat ion j'at terns of both Rayleigh and 
Love v/.r’' are two lul l'd with a nodal direction 
in E iO'* i; for R.iyleirh waves and in N 10*^ I' for 
Love v.MVi Thi.s radiation jrat t ern precludes any 
dcnible iiuple meclianism. The radiation pattern, 
the ini I ia> )>has<', the relalivc^ly larp.c' amplitude 
ratio et t.ove to Rayleii'.h waves and the existence 
of cJeai nodes of the Rayl i' igh- wav(' radiation 
pattern snpi'.est that tin' sie.u'ce is rejneseiit C'd by 
an almost hiu ixontal (Ic-ss than 15^ from tlie hor- 
izontal) siiip.le force' pointed toward S 10^' W; no 
iscrtrojiii eomponent is allowed. The surface v;ave 
spectra fall off very rapidly at ]>eriods sliiuter- 
than r.D si c .'.uitgest iirp a very slow sour ce' ju ■.> 
ccss. Althunp,!) th(' details of the source time 
history cinild not bc' det c'l niine'd , a smooth lu'il- 
shapc'd I ime function having a total dural ice of 
about 160 se'c is considered appropriate' on the 
basis of compa I' i f'.oir of synthetic and the' observed 
sci snioj’.r '""i' • ’fbe pc>ak veilue of the fore e- is 
about 10^'^ dynes. The result can be interpre'ted 
in terms of a northward lateral blast , a north- 
ward l.ind;: I iele or combination of both. The mag- 
nitude ;ni(l the' time history of the fore-e provide' 
importairt c lUi;'. t r ,i int s on the physica'. mechairlsm 
of the' e' nip t i on . 
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Seismicity and Volcanic Rock Series in Island Arcs 

Ulroo Kanamori and Marius S. Vassiliou 
Selsmologlcal Laboratory 
California Institute of Technology 
Pasadena, California 91125 

Plate subductlon Is generally considered to be the ultimate cause of large earth- 
quakes and volcanism at island arcs. The activity of large earthquakes can be 
related to various parameters of plate subductlon. For example, the magnitude, M^, 
of the largest characteristic earthquake for each subductlon zone can be expressed 
reasonably well by a linear function of plate convergence rate and the age of the 
subducting slab (Ruff and Kanamorl, 1980). On the other hand, various parameters for 
arc volcanic rocks have been associated with other characteristic features of Island 
arcs. Among them are the ratio CA/TH (Calc-Alkali Ser ies/Tholei i t Ic Series) 
(Miyashlro, 197A), Silica Index 0 (Sugimura, 1968), Alkali (K20, Na20 etc.) (Kuno, 
1966; Dickinson and Hatherton, 1967; Dickinson, 1968; Nielson and Stoiber, 1973; 
Suglsaki 1976; Ui and Aramaki, 1978) and H 2 O contents (Sakuyama, 1981). These 
parar.eters have been correlated with the depth of the Wadati-Beniof f zone; crustal 
thickness, Bouguer gravity anomalies and the plate convergence rate. We investi- 
gated various relations between seismic activity measured by My and the parameters 
for volcanic rocks to clarify the relation between the stress state at subductlon 
zones and magma genesis discussed by llycda and Kanamori (1979). We found that 
CA/TH ratio and alkali content correlate better with My than with the cmstal 

thickness. Since a large value of My can be considered to be a result of ln< leasud 
tectonic stress at nlate boundaries, tl»e above correlation suggests th.at the 
regional difference in the level of tectonic stress is at least partially rc.sponsi- 
blc for the regional variation of volcanic rock scries. 
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• SUBDUCTION ZONE SEISMICITY AND STRESS IN SLABS 


M . S . Vasst ilou 
B. H. Hap.cr 
D. L. Anderson 
A. Racfsky 

(all at: Selsmological Laboratory 252-21, 

California Institute of Technology, Pasadena 

CA 91125) 

Histograms of log N (tf earthquakes) plotted 
against depth (NOAA & PUE catalogs, > A, 
196A-1980) have the following featur .s: a 

roughly linear decay down to 250-300 km., 
followed by cither (1) a complete cessation or 
(2) a minimal level in turn followed by a 
resurgence. There Is a cessation beyond 650-700 
km, perhai>s associated with the 670 km seismic 
discontinuity. The dei’rce to which the initial 
decay is lin< ,ir varies among regions, ns docs 
the slope of the decay. Slopes correlate only 
weakly witli subduct ion parameters such ns 
convergence rate*. 

A model of a viscous slab sinking vertically 
in a less viscous medium with a rigid bottom 
boundary at 070 km yields a slab stress 

magnitude wh i . h decreases linearly with dcjith 
down to ~ 250 km, and then increases. This 
simple model, with an assumption that N varies 
exponentially as the stress, can explain the 
Important features of the seismicity curves. 
More complex models including dipping slabs and 
various boundary conditions differ in details, 
but there is a broad class of models which 
preserves. ttie essential characteristics 

mentioned above. 
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A DYNAMICALLY CONSISTENT CALCULATION OF GEOID 
ANOMALIES IN A VISCOUS PLANET WITH IMPLICATIONS 
FOR CEOID INTERPRETATION 

Ma^r k A. R i r h n v^ds 

Bradford U. llagor (Botli at: Seismological 

LfiboraLury , Caltecb, Pasadena, CA 91125) 
(Sponsor; 11. H. Hager) 

The Ifirih and other terrestrial planets are 
bertor apiaci>,i nia ted as viscous than as rigid 
spbercu;. fnlrrior density contrasts caus<:“ flow 
and stre.sst s wLicli lead to deformation of tlic 
outer suri.ii ■ and inii,rior l)Oundaries- lliese 
deformed sin ( aces cause additional geoid 
anomalic's wliieli must also bc' considered in 
evaluatin;’, (In- total g/uid anomaly resulting, 
from an int' i ior den;; i : y contrast. 

We caliul ale self-consistent geoid anomalies 
for den:;lty lontrasts in a viscous 
sel f-grav i t at ing, s|)lierc at; a function of 
spherical harmonic degree f. Free-slip boundary 
conditions ate ai'plied at the core-mantle 
boundary and no -si ip boundary conditions are 
applied at tiie surface. Additional mass 
anomalies arls.c- wlien the boundaries are deformed 
by the normal stresses. For a uniform mantle, 
defornintion of thi; boundaries leads to a larger 
contribution t f> tlie total geoid anomaly than the 
density coni i ast ilf.clf for all t, resulting in 
negative g.ioid .luonulics for positive density 
contrast r.. For intermediate f, geoid anomalies 
from density contrast.s In the lower mantle are 
an order of magpiilude larg,er than tliosc ohtaineci 
for a rig, id l.ai tli, reopening tlie question of the 
ma.ximmn dc[itii of density anomalies wtiicli can 
generate tin- observed g,eoid. 

Fot a |ilanei witli chemical or visco.si fy 
layering, tlie Ixdiavior is more coinpl ical I'd . Tlie 
sig,n and in ij.n i tude of tlie geoid anomaly depend 
u))on tlie viscosity .structure, wavelength, and 
depLli of tile density anomaly. Wc conclude' that 
i n t e r pr e t . 1 1 ! oil of g,ec>id data in terms of den.sitj 
contrasts leiniies gross assusip t i onf. about the 
mantle's rle'o 1 c)g,v and composition. 
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THE OCT. 1980 EARTHQUAKE SEQUENCE IN A SEISMIC 
GAP NEAR THE LOYALTY IS. , NEW HEBRIDES 

John Vldale (Selsroologlcal Laboratory, 

Calltornia Institute of Technology, Pasadena, 
California 91125) 

Hlroo Kanamori (Same) 

Four large earthquakes occurred in a seismic 
gap of category 2 (McCann et al. , 1979) near the 
Loyalty Is. in the New Hebrides. At 3:25 on 
Oct. 24, 1980, an event with Mg»6.7 initiated 
the sequence. Three events, Mg«6.7, 7.2, and 
6.5 followed this event on the next day. 

We investigated this sequence by using the PDE 
data, first-motion data from WWSSN records, and 
long-period surface waves recorded by IDA, SRO, 
and ASRO stations. The first-motion data 
constrained one of the nodal planes. With this 
constraint, inversion of Rayleigh and Love-wave 
spectra at 256 seconds determined the other 
nodal plane. Tlie mechanisms of all four events 
are almost pure thrust on a plane dipping 20° 
and striking parallel to the local strike of the 
New Hebrides trench. The first-day aftershocks 
indicate an Initial rupture zone of about 2,000 
km^ which is consistant With the estimated 
seismic moment of 3 x 10^ dyne-cm. During the 
next two days, the aftershock activity expanded 
to an area of 10,000 to 20,000 km in the 
directions both along and perpendicular to the 
trench. This pattern suggests that the initial 
rupture zone represents a zone of increased 
strength (i.e. asperity), and the stress change 
due to failure of this asperity subsequently 
migrated outward. During the 2 year period 
before the main event, seismicity in the initial 
rupture zone had been very low except near the 
point where the main shock sequence initiated. 

A very tight clustering of activity occurred 
there. This pattern indicates gradual stress 
concentra'.ion near the asperity which finally 
failed during the main-shock sequence. 

Since the physical condition of the fault zone 
is expected to persist in time, the present 
results would help predict the behavior of this 
gap during the next sequence. 
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AN ARRAY STUDY OF THE UPPER 
MANTLE BENEATH WESTERN MEXICO 
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Marianne C. Walck (Selamological Laboratory, 
California Institute of Technology, Pasadena, 
California 91125) 

Don L. Anderson 

The Caltech - USGS Southern California Seismic 
Network (SCARLET) is favorably located to record 
short-period P-wave arrivals from earthquakes in 
Mexico and central America (Az “ 330^15°, 

A < 40°) which have traversed upper mantle paths 
beneath tectonically active regions. The large 
number of moderate (5.0 < < 6.0) earthquakes 

occurring In this area since 1977 provide an 
excelleni data base for the investigation of the 
Gulf of California's upper mantle velocity 
structure. Record sections have been 
constructed for 29 events ac distances of 
8°-40°; each section has a station spacing of 
~10 km and extends 4.5° in distance. Arrivals 
from the triplications associated witli the two 
major upper mantle discontinuities are clearly 
visible for many events. Times of all prominant 
phases were used in direct determinations of 
dT/dA for each earthquake. Because of severe 
near-receiver structural effects, we first 
applied station corrections derived from 

relative arrival times of distant (A > 30°) 

events along the same azimuth. The resulting 
dT/dA profile has little scatter. Abundant data 
between 18°and 25° strongly constrain the 

velocity gradient between the '400 km' and 
'650 km' discontinuities; in this range the 
data is similar to that of King and Calcagnile's 
(1976) model for northern Eurasia. 
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